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In the current paper, a method is developed to diagnose potential electrical 
defects attacking doubly fed induction generator in Moroccan wind turbines 
firms. The proposed method is based on stator and rotor currents Lissajous 
curves analysis. Firstly, we focus on modeling of a non-defected wind 
conversion systems based on mathematic model created in Matlab Simulink 
which is able to reflect the behaviour of the wind turbine during 
asynchronous generator defects-free operation. After that, an indirect stator 
field vector oriented control is applied to obtain the wind system 
performance. Finally, stator and rotor currents Lissajous curves are analyzed 
in case of a non-defected generator that represents the system reference 
curves for diagnosing defects. The simulations had been realized by Matlab 
Simulink. Their results showed the effectiveness of the proposed method. 

Copyright © 2019 Institute of Advanced Engineering and Science. 

All rights reserved. 


Corresponding Author: 

Fatima El Hammouchi, 

Department of Eleetrieal System, 

Hassan 2 University, 

National Higher Sehool of Eleetrieity and Meehanie, Casablanea 8118, Moroeeo. 
Email: elhammouehi.fatima@ gmail.eom 


1. INTRODUCTION 

It was extremely pressing, both for environmental and eeonomie aspeets, to look for alternative and 
elean energies in order to satisfy growing eleetrieity needs, and reduee the C02 emission [1]. For these 
reasons, many eountries start exploring their natural and sustainable resourees sueh as wind. Aetually, 
eonsiderable wind potential is loeated in lot of parts of the world. This potential is eompletely free, 
inexhaustible and renewable. However, the wind energy eost remains high due to the expensive operational 
and maintenanee eosts (O&M) of wind turbines. Indeed, these big struetures are usually situated in remote 
loeations that are hard to aeeess. So, the maintenanee eosts ean reaeh from 10% to 30% of the wind turbines 
total ineome [2]. These eosts eould be redueed by monitoring eontinuously the system health [3], espeeially 
the generator. This erueial part of wind turbine is responsible of eleetrieity produetion and its transmission to 
the grid. 

In faet, the Doubly Fed Induetion Generator (DFIG) is the most used in the wind turbine variable 
speed industry [4]. It offers exeellent operational and eontrol features that make their integration with power 
grids easy and effeetive [5]. Mueh researeh deals with different methods to improve the eleetrieal energy 
performanee and produetion from wind [6]. Most of it develops eontrol lows for DFIG to extraet maximum 
eleetrieal power from wind turbine [7, 8]. But, during their operation, the DFIG ean be exposed to diffieult 
eonditions or manufaeturing defeets [9, 10], whieh ean immediately paralyze eleetrieity produetion and lead 
to a loss of strueture [11]. Thus, the defeets must be predieted on time. 
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A lot of works have been done to inspeet DFIG faults. Some researeh is using available teehniques 
like Supervisory Control And Data Aequisition SCADA [12]. Others use data derived from eondition and 
Struetural Health Monitoring (SHM) to deteet wind turbine faults by using deeision tree learning algorithms 
with big data [13]. 

Henee, the generator diagnosis defeets provides early deteetion of imminent faults, and minimizes 
wind turbine downtime [14]. Also, it antieipates serious damage affeeting wind energy system by predieting 
any defeets, stopping their propagation and saving time and money. In order to deteet defeets and their 
eauses, the diagnostie system has to use data from normal operation and eompare it with defeets 
eharaeteristies. So, a new method is proposed for having reliable referenee system in ease of non-defeeted 
DFIG so as to prediet the prospeetive defeets. 

In the first part of this researeh, a deseription of the essential wind energy eonversion systems 
eomponents is presented. In the seeond part, wind turbine behavior during defeets-free operation of DFIG is 
treated using mathematieal equations to ereate model in Matlab Simulink. In the last part, stator and rotor 
eurrents Lissajous eurves are developed as a method for DFIG diagnosis defeets. The simulation results, 
obtained by Matlab Simulink, show the effieieney of the proposed approaeh. 


2. WIND CONVERSION SYSTEM MODELING 

The study will be about two Moroeean wind turbine firms that eontain horizontal-axis HAWT using 
DFIG. In faet, the DFIG is an asynehronous maehine fed from the stator and rotor. The rotor of DFIG is 
eonneeted to the grid through the transformer and the voltage souree eonverters “baek-to-baek”: the Rotor 
Side Converter (RSC) and the Grid Side Converter (GSC). A eapaeitor is plaeed between two eonverters to 
reduee the voltage ripples. The RSC eontrols features of the aetive and reaetive powers, the developed torque 
and the rotation speed. So, the GSC eontrols the DC bus voltage and the grid faetor. 

2.1. Wind turbine modeling 

The wind speed used in the model is a sum of sinusoidal signals with different frequeneies and 
amplitudes as shown in the Figure 1 below. 


vmd speed 



Figure 1. Wind speed model 


By applying the theory of momentum and Bernoulli’s theorem, ineident theoretieal power due to 
wind is expressed by the equation below: 



Pwind • Wind power; S : Air Density; A: Sweep area swept by the blade; v : Wind speed. The Power eaptured 
by the blade is determined by the equation below: 

( 2 ) 

With 

C/A,fi) = ki.{k2.K-k,.j3-k,y‘^^ +k,.A 
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Cp(A,,P)\ Coefficient of performance 
^: Tip-speed ratio 

P : Pitch angle (=0 for maximum ) 

Wturb: Angular speed of the blade 

In order to have y^=0 for a maximum C^ : kl= 0.5872, k2=l 16, k3=0.4, k4=5,k5= 21, k6= 0.0085, 
k7=0.08, k8=0.035 [15]. The aerodynamic torque on the main shaft (slow axis) can be expressed by the 
equation: 




«5) 


7 = 



+7. 


(7) 


Tturb: Blade torque 

J: Inertia of the system brought back on the turbine fast axis 
Jturb: Turbine inertia 
G: Gearbox ratio 
Jg: Generator inertia 

The fundamental equation of dynamics is applied to determine the evolution of the DFIG 
mechanical speed which is presented in equation below: 

J^^ = T,-Te«-fy^n,ec («) 


Wmec: DFIG mechanical speed 
Tg: Aerodynamic torque on the the turbine fast axis 
Tern* Electromagnetic torque 
f: Friction coefficient 

The equations presented above, are used to create the wind turbine Matlab-Simulink. We suppose 
that the frictions are neglected. The system performance is showed in Figures 3, 4, 5 and 6. 

2.2. Mdeling of doubly-fed induction machine 

The asynchronous machine equations for the stator Us and the rotor Ur voltages in Park reference 
frame (d, q) are [16]: 


d (/> as 

dr = + Rr-i dr “ ® s " ® r ) 

d ^ qr 

U qr = ) 


(9) 


With: cOj. -co^ -P.w (10) 

Ug (d,q), Ur (d,q): Stator and rotor voltage Park components respectively; 

(d, q), 4 (d, q): Stator and rotor flux Park components respectively; 
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ig (d,q), ir (d,q): Stator and rotor currents in Park components respectively; 

Rs, Rr: Stator and rotor resistanees respeetively. The flux equations ean be written as follows: 


= L, .i* + M .irf, 

= L, + M 

^ dr = -^dr + ^ ds ( 1 ^) 

= Lr J qr + ^ qs 

L ^ = I ^ - M g , 

Ls, Lr : Stator and rotor eyelie induetanees respeetively; 

Is, lr : Stator and rotor leakage induetanees respeetively; 

Ms, Mr : Mutual induetanees between stator phases and rotor phases respeetively; 

M : Maximum mutual induetanee between stator and rotor stages; 
p represents the DFIG pole pairs number. The eleetromagnetie torque expression is: 

'^em ~ P( ^ds "^qs ~ ^qs "^ds ) ( 1 ^) 

2.3. Field oriented control 

The purpose of orienting the stator field is to perform the DFIG aetive and reaetive powers by 
eontrolling dynamieally and separately the flux and the torque. That means orienting stator flux along the 
direet axis in order to obtain: = 0 and ^g = <pgd . From (11) and (12) we get: 


Ls 

. . 

- —-v 

T A ■ ^ A ■ 

Tern =P9s-^qs 


(13) 


In general, the stator resistanee is negleeted for large power equipment used in wind turbine. It is 
supposed that the grid is stable and the field is eonstant. The stator voltage park eomponents (9) beeome: 




The DFIG aetive and reaetive powers ean be written as follows: 
J ^ ds ’^ds ^ qs ’^qs 

[Ss “ ^qs ’^ds ~ ^ds ’^qs 


(15) 


By using (14), the aetive and reaetive powers expressions beeome: 




Ls 


Qs = 




M 


Ls 


(16) 


So the aetive power is independent of the reaetive power. It is eontrolled by the quadrature rotor eurrent iqr. 
The reaetive power is eontrolled by direet rotor eurrent idr. 
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=(7?, +S(4 ——))4 - -^)-V 

Ls Ls 

M 

u,. = (K + siL^ - —)K + (4 - —K 

There are two approaches of the vector control: direct and indirect. The direct one is simple to 
implement [17] but not the most efficient [18, 19] because it ignores the terms of coupling. The indirect 
method is being more commonly used [20] because it allows operating more easily throughout all the speed 
range [21]. 

In order to have performing system, the wind turbine model is implemented in Matlab Simulink by 
using the indirect method (IFOC). The orientation of the DFIG stator flux is chosen according to the direct 
axis. The overall pattern of the wind system model is presented in Figures 7, 8, 9 and 10. 


3. DIAGNOSIS OF NON DEFECTED WIND CONVERSION SYSTEM 

Unfortunately, wind turbine faults occur once it starts working. Indeed, Wind energy has its 
congenital defects due for example to the harsh operating environment and the highly turbulent wind speed 
[22]. But, under whatever operation conditions, the stability and reliability of wind turbine are requisited to 
the grid. Hence, the diagnosis of wind turbine should be realized. In this part, we diagnose a non-defected 
system to have a reliable reference. 

3.1. Review of diagnosis defects essential methods 

For minimizing downtime, predicting in time imminent failures and improving maintenance 
schedule, several techniques are applied for monitoring continuously the wind energy conversion system 
health [23]. Indeed, online diagnosis defects allows an early detection of the structure health degradation and 
prevents incipient faults in the wind energy conversion system which leads to save both money and time. 
Multitude of internal and external defects attacks frequently wind turbines once its start working [24] 
especially the generator. Indeed, DFIG undergoes some common defects such as generator excessive 
vibration, generator overheating, bearing overheating, abnormal noises and insulation damage, etc [25]. 

Much research on generator faults deals with methods of process simulation and algorithm 
construction which are used in condition monitoring to detect and anticipate defects [26]. Other techniques 
use time-frequency analysis tools, such as wavelets, to derive a fault detection signal [27]. Some methods are 
also exploring data provided by the supervisory system and applying many data-mining algorithms to 
develop models for predicting occurrence of possible defects [28]. These different techniques are widely used 
in control system for reducing the wind turbine downtime and minimizing operational and maintenance costs. 

3.2. Method developed for DFIG diagnosis defects using Lissajous curves 

The Lissajous curves are studied first by the American mathematician Nathaniel Bowditch in 1815. 
The curves, in 1857-1858, were investigated independently by the French mathematician Jules-Antoine 
Lissajous. The French mathematician used a narrow stream of sand pouring from the base of a compound 
pendulum to produce the curves [29, 30]. 

Lissajou figures are created by the combination of two sine waves: the x signal is described by one 
sine wave, and the y signal is described by another sine wave. So by expressing y as a function of x, 
interesting patterns are created. Recent developments in the field of Lissajous curves are realized in different 
domains such as the rotating machinery field [31], and faults diagnosis [32]. 

Lissajous curves are described by the parametric equations below [33]: 

X = A sin(at + a) y = B sin(bt) (18) 

If the two waves have exactly the same frequency (a=b), the same amplitude (A=B), with an angle a=n/2, the 
parametric as shown in 18 become: 

X = i4cos (at) y = Asin(at) (19) 

x^ + = A^ ( 20 ) 

Then the resultant output will be a circle. In this view, Lissajous curves are applied in this paper to the DFIG 
stator and rotor three-phase currents. For this purpose, it is necessary to use the Concordia transformation 
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[34] in order to obtain the stator and rotor eurrents Coneordia eomponents (isa, ira) and (isp, irp) from the 
three-phase stator basie eurrents isabc and the three-phase rotor basie sueh as: 


ia 

- H 

.W. 

3 


■ 1/2 

- 1/2 


ia 



* 

ib 

2 

2 . 


Ac. 


( 21 ) 


Matlab Simulink software is used for simulating a non-defeeted system. The results are presented in Figure 
13 and Figure 14. These eurves are eonsidered as a referenee to whieh the defeeted system eurves will be 
eompared in future works. 


4. RESULTS AND ANALYSIS 

The results below represent the simulation of the wind turbine model in Matlab Simulink at the first 
thirty seeonds. The results are obtained by ehoosing the following values on Table 1. 


Table 1. Model parameters 


Turbine 

parameters 

Coeffieients for 
maximum Cp 

Maehine 

Parameters 

S =1.225 kg/m^ 

kl= 0.5872 

Rs = 0.012 n 

A = 3903.625 

k2=116 

Rr = 0.021 n 

R = 35.25 m 

k3=0.4 

Ls = 0.0137 H 

G=90; 

k4=5 

Lr = 0.0136 H 

j=1000; 

k5=21 

M = 0.0135 H 

f=0.0024; 

p=o° 

k6= 0.0085 
k7=0.08 
k8=0.035 

g =0.003 


The wind system is exposed to variable wind speed from 8 (m/s) to 14 (m/s) as shown in Figure 2. 
From (1) we obtain Figure 3 whieh represents the wind power Pwind- Then, by using (2) we get the power 
eaptured by the blade Pturb as shown in Figure 4. The Figures 5, 6 and 7 show the system performanees. In 
deed, we ean see that the tip speed ratio is almost eonstant as illustrated in Figure 5. The eoeffieient of 
performanee reaehes its maximum value whieh is 0.57. 



Figure 2. Wind speed (m/s) 


Pwind (w) 




25 Time(s) 30 




Figure 3. Pwind 



Figure 4. Pturb power eaptured by the blade 



Tijnr (t) 


Figure 5. Tip-speed ratio 
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Figure 6. Coefficient of performance 


Figure 7. Mechanical speed of DFIG 


The stator flux of DFIG is oriented according to the direct axis that is why the machine stator 
indirect flux, in Figure 8, reaches zero and the stator direct flux is constant as shown in Figure 9. In addition, 
the indirect field oriented control is used to control active and reactive powers. A reference for active power 
Psref is chosen to have a form of steps (3 steps). The Figure 10 shows a comparison between power reference 
and real active power. It is clear that the real active power follows the power reference. The system real 
reactive power is also around the reactive power reference which is zero Figure 11. Stator three phase 
currents are presented in Figure 12 from which shows that DFIG three-phase stator currents isabc have a 
perfect sinusoidal shape. By using (15), the stator currents are influenced by Psref For this reason, three 
variation zones of stator currents are observed. These figures show that the wind system model implemented 
in the Matlab Simulink react rapidly and has a satisfactory dynamic. 



Figure 8. Indirect stator field of DFIG 


Figure 9. Direct stator field of DFIG 



• * ^ '* ** ■’' riiw(*} " *g ^ A Twwti) * 

Figure 10.Comparison between Psref and Ps Figure 11. Comparison between Qsref and Qs 



Figure 12. Stator currents phases a, b and c 


The Figures 13 and 14 expose Fissajous curves of a DFIG sator and rotor currents during free- 
defects operation. Each curve, has three circular shapes passing from the small circle to the biggest one. 
Theses three circles are obtained because the Psref is chosen as a signal with three steps. In deed, at the 
simulation starting time and after application of the active power reference signal first step Psref, the first 
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circle is obtained. At 0.3 second, another step is added to the first step. Thus, the second circle is drown. 
Finally, at one second, we apply the third step that lead us to get the third circle. 

The rotor currents Lissajous curves are showed in Figure 13. It has a 6 petals flower shape. This 
shape is due to the DFIG rotor supply by the inverter with 6 switches. The DFIG currents Lissajous curves 
will give an idea, in fure works, about the kind of the defects attacking wind conversion system. So, when the 
defects occur, the Lissajous curves will have other shapes. 



Figure 13. Rotor Lissajous curves 


Figure 14. Stator Lissajous curves 


5. CONCLUSION 

In the first part of this paper, the wind conversion systems modelling based on mathematical 
equations is presented. Then, an indirect vector control stator field oriented is applied to increase the 
performance of the wind conversion chain. The system is simulated by Matlab Simulink during defects-free 
operation. In the second part, a method was developed to predict the potential electrical defects attacking 
wind turbine using Lissajous curves of the DFIG stator and rotor currents. 

The obtained results are considered as a reliable reference system for DFIG diagnosis of Morocan 
wind turbines parks. To this reference system, the wind conversion energy behaviour will be compared in the 
case of the defects. The next work will be devoted to modelling and simulation of defected wind energy 
conversion system. Lissajous curves and frequency spectrum will be used in order to analyse and diagnose 
different defects attacking the two Moroccan wind turbines parks. 
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